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Abstract 



s^ , We revisit leptogenesis in the minimal non-supersymmetric type I see-saw mech- 

CN \ anism with two right-handed (RH) neutrinos, including flavour effects and allowing 

f^ \ both RH neutrinos A'^i and A'^2 to contribute, rather than just the lightest RH neu- 

^O \ trino A^^i that has hitherto been considered. By performing scans over parameter 



space in terms of the single complex angle z of the orthogonal matrix R, for a range 
of PMNS parameters, we find that in regions around z ~ ±'7r/2, for the case of a 
normal mass hierarchy, the N2 contribution can dominate the contribution to lepto- 
genesis, allowing the lightest RH neutrino mass to be decreased by about an order 
of magnitude in these regions, down to Mi ~ 7 x 10^*^ GeV for initial thermal N2- 
abundance, with the numerical results supported by analytic estimates. We show 
that the regions around z ~ ib7r/2 correspond to light sequential dominance, so the 
new results in this paper may be relevant to unified model building. 



1 Introduction 

Current low energy neutrino data with bi-large mixing can be minimally explained within 
the non-supersymmetric type I seesaw mechanism [l| with only two RH neutrinos [2|. 
This can be regarded as a limiting case of three RH neutrinos where one of the RH 
neutrinos decouples from the see-saw mechanism either because it is very heavy or because 
its Yukawa couplings are very weak. With only two RH neutrinos it is straightforward 
to see that the lightest left-handed neutrino mass has to vanish. Since the number of 
parameters is greatly reduced, this model has attracted attention in connection with the 
possibility to test it, especially when successful leptogenesis |3| is required in addition 
to the constraints from low energy neutrino data [^. The number of parameters (11) is 
however still high enough that the parameter space cannot be yet over-constrained by 
combining together low energy neutrino experiments and successful leptogenesis. Due to 
its minimality, there has been a great deal of attention paid to the two RH neutrino type 
I non-supersymmetric see-saw mechanism, (i) in the unflavoured case, then (ii) including 
flavour dependent effects, as follows. 

(i) In the unflavoured approximation, in p] it was shown that, imposing double texture 
zeros in the neutrino Dirac mass matrix, it is possible to make connections between the 
sign of the baryon asymmetry of the Universe and the sign of CP violation in neutrino 
mixing. Similar results also apply when the two RH neutrino model is regarded as a lim- 
iting case of three RH neutrinos when the heaviest RH neutrino of mass M-i ^10^"^ GeV 
decoupler from the see-saw mechanism g 3. In B a systematic leptogenesis study ot 
this model was performed, also with texture zeros in the neutrino Dirac mass matrix, con- 
sidering lepton flavour violating processes within supersymmetric scenarios. Leptogenesis 
with two RH neutrinos has been also studied beyond the hierarchical limit, obtaining 
the precise__conditions both for the degenerate limit and for the hierarchical limit to be 
recovered 



(ii) A first analysis that included flavour effects [10|] was presented in [ll|. The close 
connection between the CP violation for leptogenesis and the observable leptonic Dirac 
CP phase in models with two texture zeros and in the limit of a nearly decoupled RH 



neutrino, which corresponds to the two RH neutrino case, has been studied in 121]. In 



13l | a flavoured analysis of leptogenesis within the two RH neutrino model has shown that 
for an inverted hierarchical spectrum and if the condition —sin ^13 cos 5 > 0.15 applies, 
then the model can work only in presence of CP violating Majorana and Dirac phases. 
In [iJ] the two RH neutrino limit was considered within a leptogenesis scenario when CP 
violation is uniquely stemming from the Dirac phase. The leptogenesis lower bound on 



Ml in the two RH neutrino model has been studied in 16|, showing that in the absence 



of PMNS phases and for inverted hierarchy the bound is much more strin gent , confirming 
the importance of the PMNS phases for this case, as also pointed out in |l3| . 

In this paper we go beyond the above analyses of leptogenesis in the minimal non- 
supersymmetric type I see-saw mechanism with two RH neutrinos, allowing both RH 
neutrinos A'^i and A'^2 to contribute, rather than just the lightest RH neutrino A'^i that 
has hitherto been considered. We also include full flavour effects which are crucial for 



A^2 leptogenesis 17| and emphasise the role of a self-energy contribution to the flavoured 
CP asymmetries which is frequently ignored for A'^i leptogenesis in the hierarchical limit, 
but which will prove to be crucial for successful A'^2 leptogenesis in the hierarchical limit 
assumed in this paper. By performing scans over parameter space in terms of the sin- 
gle complex angle z of the orthogonal matrix R, for a range of PMNS parameters, we 
find that in regions around z ~ ±7r/2, for the case of a normal mass hierarchy, the 
A^2 contribution can dominate the contribution to leptogenesis, allowing the lightest RH 
neutrino mass to be decreased by about an order of magnitude in these regions, down 
to Ml ~ 7(10) X 10^*^ GeV in the case of initial thermal (vanishing) A'^2-abundance. In- 
terestingly, these regions correspond to so-called light sequential dominance, in which 
Ni dominantly contributes to the atmospheric neutrino mass nis, while N2 dominantly 
contributes to the solar neutrino mass m2. 

In Section 2 we set up the general notation for leptogenesis in the two RH neutrino 
model. In section 3 we solve the Boltzmann equations finding an analytical solution for 
the final asymmetry. In section 4 we express the CP asymmetries within the orthogonal 
parametrization. In Section 5 we show the allowed regions contrasting the A'^i production 
and the N2 production and showing that new regions appear thanks to the A'^2 contribution 
to the final asymmetry. In Section 6 we show that these regions correspond to light 
sequential dominance. In Section 7 we draw the conclusions. 

2 General Set up and notation 

The Yukawa part of the Lagrangian in a SM extension to include heavy RH neutrinos is 
given by, 

-Cy = YJHIr + Y.LHNr + h^^MNn + h.c, (1) 

where L and H are the left-handed lepton doublet and Higgs doublet respectively, Ir the 
RH charged singlet and Nr the RH neutral singlet. Y^. and Y^ are the Yukawa couplings 
and M the RH Majorana neutrino mass matrix. In the above equation H = —ia2H*. 



After electroweak symmetry breaking we get the Dirac mass matrix m^, = Y^v, where v is 
the vacuum expectation value of the Higgs doublet. If we consider n generations of heavy 
RH neutrinos N^, then the Dirac mass matrix tti/j is a 3 x n matrix and the Majorana 
mass matrix M is a n x n matrix. The (3 + ri) x (3 + n) neutrino mass matrix turns out 
to be, 

mD\ ( ^l 



In the see-saw limit, M ^ m^), the spectrum of neutrino mass eigenstates splits in two 
sets: n very heavy neutrinos, Ni, . . . , N^ respectively with masses Mi < M2 < . . . < M„ 
and almost coinciding with the eigenvalues of M, and 3 light neutrinos with masses 
nil < ""^2 < ""^3 for normal hierarchy (NH) and 7723 < rrii < 7712 for inverted hierarchy 
(IH), the eigenvalues of the light neutrino mass matrix. 

Once the n heavy RH neutrino fields get integrated out from the theory, one obtains 
the 3x3 light neutrino mass matrix, up to an irrelevant overall sign, as 

m^ c::i mo M~^ m^, (3) 

where we neglected terms higher than 0{M^^). The heavy neutrino mass matrix is 
approximately given by M. In this paper we assume M3 ^ 10^^ GeV so that N3 decouples 
from the seesaw and effectively a two RH neutrino model {n = 2) is recovered. 

Let us now introduce the relevant general quantities for leptogenesis. First of all notice 
that in addition to M3 ^ 10^^ GeV, we will also impose the condition M2 ^ 10^^ GeV. In 
this case lepton flavour effects have to be taken into account at the asymmetry production 
from A^i and N2 decays 10|. 

This can be done calculating the total B — L asymmetry as the sum, 

Nb-l= Y. ^^^^ (4) 

a=e,iJ,,T 

of the flavoured asymmetries defined as Aq = B /?> — La. Notice that Nx indicates any 
particle number or asymmetry X calculated in a portion of co-moving volume containing 
one heavy neutrino in ultra-relativistic thermal equilibrium, i.e. such that N'^_[T ^ Mj) = 
1. The baryon-to-photon number ratio at recombination is then given by 

^B = a.ph^=.0.01iVi_^, (5) 

7 



to be compared with the value measured from the CMB anisotropics observations 221] 

r7g^^ = (6.2±0.15)x 10"^°. (6) 

4 



The Ni RH neutrinos' decay widths are given by 

r, = r, + r, = ^ — -^^ , 7 

where Fj and Fj are respectively the Ni decay rates into leptons and anti-leptons (at zero 
temperature). The key quantities encoding the main features of the kinetic evolution are 



the decay parameters defined as [l8|, [l9| 



j^ Ti nii ^ {m'^mD] 

Kj = --—:z r-TT = — , where m,- ~ 



H{T = Mi) m. ' ^ M, 



are the effective neutrino masses 20| and ?n^ is the (SM) equilibrium neutrino mass defined 



by [19|,|21| 



16 vr^/'^/^ 



V 



2 



^ ZZ:: ^^£i!i__ ~ 1.08 X 10~^eV. (9) 

3v^ Mpi ^ ' 

The flavour composition of the lepton flavour quantum states produced by the Ni decays 

can be written as, 

Ki) = 5^ C,„ 14) , C,^ = {L\^i) (^ = l,2,a = e,/i,r) (10) 

a 

and 

\Ii) = J2 ^ic, 14) , Ci^ = {ia\i^) , it = l,2,a = e,fi,T), (11) 

a 

where we notice that in general the flnal anti-lepton states produced by the N2 decays 
are not in general the CP conjugated of the flnal lepton states and therefore, in general. 



Cia 7^ Cia [1Q|. Introducing the branching ratios Pia = |CjQ,p and Pjq, = |CjqP respectively, 
i.e. the probabilities that a lepton or an anti-lepton is measured in the a light lepton 
flavour eigenstate, we can deflne the flavoured decay rates Fj^ = Pia Fj and Fjq = Pia Fj. 
In a three flavoured regime, when the produced lepton quantum states rapidly collapse into 
an incoherent mixture of flavour eigenstates, the Fj^ can be identified with the flavoured 
decay rates into a leptons, F(A^j — )■ (j)^ la) , and anti-leptons, T{Ni — )■ ^L), respectively 
and the Pia and Pia with their branching ratios. The total flavoured decay rates (at zero 
temperature) are then given by 

^ la — ^ la 1^ ■'- la o 9 ' \ J 



Correspondingly, the flavoured decay parameters are deflned as 



"'-H{T = Mi) Mim, -"" '' ^ ' 
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where Pj^^ = {Pia + Pia)/'^ are the tree level branching ratios. 
We can then define the fiavored CP asymmetries as 



^iry — 



r - r 



SO that for the total CP asymmetries one has 

r, - f , 



Si = 



r. + Ti 



E' 



(14) 



(15) 



In the fiavour basis, where the charged lepton Yukawa matrix Yg and the RH neutrino 
mass matrix M are both diagonal with real and positive eigenvalues, the neutrino Yukawa 
matrix Y^, will have elements denoted as hai- The fiavored CP asymmetries can then be 
calculated using [23 1 



lQii{Wh)i 



Y^U^K^KAh^h).,] 



j¥=i 



^(Xj/Xj) ^ 2 

^/x~/x'i 3(xj/xi - 1) 



—lm[hl^haj{h^h)ji] 



where we defined x,- 






and 



2 



(1 + x) In 



1 + a; 



X 



X 



1 — X 



(16) 



(17) 



The second term on the right-hand side in the expression for Eia due to the self-energy 
diagram has so far been neglected in studies of A'^i leptogenesis within the 2 RH neutrino 
model except in 16| where however the 2 RH neutrino model was not the main focus. 
In [ij] it was noticed that this term could play a relevant role in the calculation of the 
heavier A'^2 RH neutrino asymmetry, and it will be included in our analysis l • 
If we also introduce the variables 



Ml 

— -- and 
T 



Xi z , 



:i^ 



the decay and the washout terms can be written respectively as 



DAz) 



KiXiZ ( — 
>7i 



and Wi{z) = jKi y/xllCi^Zi) z\ 



3 

i 1 



(19) 



where the averaged dilution factors can be expressed in terms of the Bessel functions, 

(1/7.) = lC^{z{)/lC2{zi). 



^It was discussed recently in 15| that this term is generically dominant in scenarios where two RH 
neutrinos form a quasi-Dirac pair. Its size is related to the non-unitarity of the Icptonic mixing matrix, 
caused by an effective dimension 6 operator. Since this scenario implies Mi ~ A/2, it will not be studied 
here. 



3 An analytical solution for the final asymmetry from 
Boltzmann equations 

After having set up the notation and framework, in this Section we write down the 
Boltzmann equations and give an analytical solution for the final asymmetry. We will 
impose throughout the paper M2 > 3 Mi, so that the hierarchical (non- resonant) limit 
holds and ^ — 1 (cf. eq flT7|) ) is a good approximation [9|. In this limit the production and 
the wash-out from the A'^2's and the production and the wash-out from the A'^i's can be 
treated as two separate stages. In a first stage, the asymmetry is produced and washed- 
out by the A'^2 processes, while in a second stage when the asymmetry is produced and 
washed-out by Ni processes. 

3.1 Production of the asymmetry from N2 processes 

Recall that in the SM, if leptogenesis occurs at temperatures T ~ Mi, where Mi is 
the mass of the lightest RH neutrino, then one has to distinguish two possible cases. If 
10^ GeV <C Ml ^ 10^ GeV, then charged /i and r Yukawa interactions are in thermal 
equilibrium and all flavours in the Boltzmann equations are to be treated separately. 
For 10^ GeV ^ Mi ^ 10^^ GeV, only the r Yukawa interactions are in equilibrium 
and are treated separately in the Boltzmann equations, while the e and /x flavours are 
indistinguishable . 

In the case of A'^i leptogenesis, it is well known that the dominant contribution from 
the first term on the right-hand side of eq. (TT6|) to the flavoured CP asymmetries Sia in 
the eq. (TT6|) is bounded from above, leading to a lower bound Mi > 10^ GeV in order 
for the CP asymmetries to be sufficiently large. We will see later that this conclusion 
also applies when the production from the next-to-lightest RH neutrino A'^2 is taken into 
account and when both terms in the eq. (fT6l) . both for ei^ and e2ai are taken into account 
as well. 

Assuming Mi > 10^ GeV, we are always in the two-light flavour regime, where only 
the tauon charged lepton interactions are fast enough to break the coherent evolution of 
the final leptons. This implies that we have to track separately the asymmetry in the 
tauon fiavour and the asymmetry in the fiavour 72 defined as the coherent superposition 
of the electron and muon components in the lepton quantum states \i2) produced by N2 
decays, explicitly 



\C2e?+\C2,?^ ''' \C2e?+\C2,\ 



1^.) = ^^^^V^ 14) + ^^^^V^ 1^.) , (20) 



and in the anti-lepton quantum states \i^^) analogously defined. 

We neglect the coupling between the dynamics of distinct flavours a 7^ /3, in the specific 
case of the two flavours r and 72 [l7j. Therefore, in the stage where the asymmetry is 
produced by A'^2 decays, the relevant kinetic equations can be written as 

dNN2 



dz 



-D2 {N, 



N2 



N' 



N2) ' 



72 



dz 

dNA^ 



dz 



-62, D2 {Nm, - iV^^) - P°^ W2 iV^, 

-e2r D2 {Nn2 - iV^l) - Pi W2 Na^ 



(21) 
(22) 
(23) 



where we defined £27 = £2e + £2fi and P27 — ^l + -^2% This set of classical Boltzmann 
equations neglects different effects that have been studied in the last years such as thermal 



masses 



2J], decoherence 25|, quantum kinetic effects 26J, momentum dependence 27f |. 



In the strong wash-out regime (i^2e + -^2^4? -^2r ^ 5) these effects give at most 0{1) factor 
corrections. 

The two asymmetries A^Aq (q^ = t, 72) freeze out at T = M2/zBa where [9| 



2.5 



0(1 ^ 10) . (24) 

At the end of the N2 production stage, at 2; ~ zb2 = max[zB{K2'y), zb{K2t)], one has 



N- 



ZB2 



-£27 i^{,K< 



27 J' 



and iV^f ^ -£2r f^{K2r) . 



(25) 



In the case of an initial thermal abundance, the efficiency factors are approximately given 
by 0,0,1281 



K{Ki, 



1 -exp ( --KiaZ^{Ki, 



(26) 



KicZji^Ki^ 

In the case of vanishing initial abundances, the efficiency factors are the sum of two 
different contributions, a negative and a positive one, explicitly 



kI = kI{K,,PI) + kI{K,,PI). 



(27) 



The negative contribution arises from a first stage where Nj^. < A^^, for zi < z^^, and is 
given approximately by 



>^'-{K..P-o 



2i 3 TT K^^ 



e 2 



- A^JV; (^cq) 



(2^ 



The positive contribution arises from a second stage where A^tv. > A^^, for Zi > z'^'^, and 
is approximately given by 



4(J^„i^°J^ 



ZB{Kia)Kio 



l-e" 



l^ia ^sC-ff ia) ^JV. (zoq) 



(29) 
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The Ni abundance at z^^ is well approximated by the expression 



eq. , N{K, 



NnA^"")^- ^^^^, (30) 



that interpolates between the limit i^j :^ 1, where z^'^ <^ 1 and N^^{z'^'^) = 1, and the 
limit Ki < 1, where ^f > 1 and iV^v^l^D = N{Ki) = 37rKi/A. 

3.2 Production and wash-out of the asymmetry from A^i pro- 
cesses 

When inverse A'^i processes start to be active at ^ ~ 1, they break the coherent evolution 
of the l^^j) quantum states [29|. We describe this decoherent effect in terms of a simple 
collapse of the quantum state neglecting decoherence effects that would be described by 
a density matrix approach. This is justified since, thanks to the condition of hierarchical 
masses M2 > 3 Mi, the A^2-decays are already switched off at this stage and they do not 
interfere with A'^i inverse processes. The stage of decoherence can then be regarded as a 
transient stage with no relevant consequences on the final asymmetry. 

Therefore, at T ~ Mi, the {i^^) quantum states can be described as an incoherent 
mixture of a £^^ component and of a i ± component. Both components are a coherent 
superposition of electron and muon flavour eigenstates. The first has a flavour composition 
given by the projection of \ii) on the e — /i plane, while the second is the projection of 
the \ii) orthogonal component on the e — /i plane. 

Let us now define the probability pu = | (^71 1^72) P- This can be calculated from the 
Dirac mass matrix as 

Within the adopted kinetic description, only the component pu {i-yj^) of \i^^) interacts with 
the Higgs in an inverse process producing A'^i. The component (1 — Pi2)|^7-l), orthogonal 
to 1^71), is untouched. Analogous considerations hold for the anti-lepton quantum states. 
In this way only the asymmetry in the flavour 71, that we indicate with A'^^ is washed 
out, while the asymmetry N^ , is not changed by Ni inverse processes. 

Therefore, under the action of A'^i decays and inverse processes, the {i^^) quantum 
states collapse into an incoherent mixture of |£^j) and \i^±) quantum states and analo- 
gously the l^ya)- Correspondingly, one has to calculate separately the two contributions 
TVa^ and Na j_ to the asymmetry, in addition to the N/^^ asymmetry in the tauon flavour. 



Pl2 = ^-^ \ ( t ^- ^^^' 



Therefore, in this stage the set of Bohzmann equations is given by 

dNN, 



dz 
dNA 

dz 

dNA ^ 

]T 

dz 
dNA^ 



- -D,{N^,-N'^l), (32) 

'' -- e^,D,{NN^-N';,\)-P^^W^NA,^, (33) 

= 0, (34) 

= £i.A(iV^, -iV^'J-AlW^i^A., (35) 



dz 

implying that Na , remains constant. Notice that effectively we have in the end, because 
of the heavy flavour interplay, a three flavour regime where the final B — L asymmetry 
can be calculated as the sum of three contributions, 

N'b-l = Nl + Ni + iVi ^ , (36) 

where 

Ni^, - Pi2e2^ti{K2T)e-'-^^''< +ei^K{Ki^), (37) 

Ni^ - {1 - Pu) e2^ k{K2^) , (38) 

Ni^ ^ e2rK{K2r)e-^'''^+eir^{Kir). (39) 



It is useful for our discussion to split the final asymmetry into a contribution from A'^i 
decays and into a contribution from N2 decays, 

Nh-L = <i^L + Nf\ , (40) 

where 

<i^^ ~ ei^ k{K^^) + Sir t^{Kir) (41) 

and 

<-L - Pl2^27 <K2-,) e-"^^^- + (1 - P12) 62^ k{K2-,) + £2r ^{K2r) C""^^^- . (42) 

In this way we clearly distinguish the effect of taking into account the asymmetry produced 
from the next-to-lightest RH neutrinos N21 which has been neglected in previous analyses. 

4 Combining the low energy neutrino data with the 
orthogonal parametrization 

In this section we recast our expression for the final asymmetry in the orthogonal parametriza- 
tion, which provides a convenient way to connect the constraints from leptogenesis to the 

10 



information from current low energy neutrino experiments and the additional parameters 
from the RH neutrino sector. 

4.1 Orthogonal parametrization for the two RH neutrino model 

The light and heavy neutrino mass matrices can be diagonalized by unitary matrices U and 
Um, respectively. Hence we have the relations U^'m^U* = Dk and Uj^MUl^i = Dm, where 
Dk = diag(mi, m2, m^) and Dm = diag(Mi, M2, M3) are diagonal matrices containing the 
light and heavy neutrino mass eigenvalues for three RH neutrinos. In the basis where Y^ 
is diagonal we identify U as the PMNS matrix. From above, one obtains, 

U^ruDM-^mlU* = Dk . (43) 

Substituting U\,jMUm = Dm in the above equation we get, 

f/t mj, U*M D-J Ul ml U* = D^ . (44) 

The R matrix is defined as |30|] g 

R = D-J^UlmlU*D^^, (45) 

where i? is a complex orthogonal matrix R^R = I. Eq. ( l45l) parametrizes the freedom in 
the Dirac matrix mD, for fixed values of U, D^ and Dm, in terms of a complex orthogonal 
matrix R. 

From eq. (I45l) . in the basis where M and Yg are diagonal, rriD parameterizes as: 

^dD-1, = UD^,R^, (46) 



M 



where D^ = diag(??7,^ , m2 , m^ ) and -D ij = diag(M^ , Mg , M3 ) for three RH 



neutrinos. To be completely explicit we can write the Dirac matrix mj;) as rriDai where 
a = e,fi,T labels the rows and i = 1,2,3 labels the columns corresponding to the three 
RH neutrinos and then expand eq. ( 146|) as: 



(A#"l/2 ;\/r-l/2 A#-l/2x ,TT 1/2 TT 1/2 tt 1/2 x 

niDeiM^ mDe2M2 moe^M^ \ Ueim-^ Ue2m^ Uesm^' \ 

niD^^M-^'^ mD^2M-^'^ nio^^M^^'^ = U^im'/^ U^2ml^^ U^^mlJ^ i?^.(47) 

mDriM^^''^ mDT2M2^''^ mDrzMl^''^ ) \Ur\vax Ur2fnl''^ Ursm^s / 

The eq. P7|) enables the Dirac matrix to be determined in terms of the completely free pa- 
rameters of the complex orthogonal matrix R, for a fixed physical parameter set U, rrii, Mj. 
For example we can scan over the parameters of R for a fixed U, mi,Mi. 



^In terms of PMNS mixing matrix R = D J—U\^mlUiUpj^[^gD J- where Upmns = U}U. 

11 



As remarked, the two RH neutrino model can be regarded as a limiting case of three 
RH neutrinos where one of the RH neutrinos decouples from the see-saw mechanism either 
because it is very heavy or because its Yukawa couplings are very weak [2|]. In our case 
we shall consider the former situation M3 — )■ 00. Then we are left with only two non-zero 
physical neutrino masses which can be identified as either m2.,m^ with mi — )■ for a 
normal hierarchy (NH), or 7711,7712 with 7723 -> for an inverted hierarchy (IH). 

For the case of two RH neutrinos of mass Mi , M2 , and two physical neutrino masses 
7712, rris, for the case of a normal hierarchy, with mi — )■ 0, 

/mz)eiM~^/^ mDe2M2^^'^\ /Ue2'ml^^ Ue^rnj'^X 
7nD^.iM-"^ 7nD^2M2^'^ = U^2my^ U^sml^'' R^ . (48) 

\771DtiMi^''^ 7nDT2M2^'^ J \ f/r2"^2 t^r3"^3 / 

For the case of two RH neutrinos of mass Mi,M2, and two physical neutrino masses 
7711,7712, for the case of an inverted hierarchy, with 777,3 ~^ 0, 

TTlDelMi^^'^ mDe2M2^^'^ \ / f/el"^}^^ f^e2"^2^^ \ 

niD^iM;'^' mD,2M2^'^ = U.^m]'^ f/^2mf R^ . (49) 

.771DtiMi^''^ 771Dt2M2^''^ / \Urim^i f/r2"^2 / 

In each case the 3x2 Dirac mass matrix rra^i^j is parametrized in terms of a 2x2 complex 
R-matrix which can be written as: 

/ cos z C sin 2; \ , ^ 

R=[ - , (50) 

\ — sm 2; (^ cos z ) 

where ^ is a complex angle and C = ±1 accounts for the possibility of two different choices 
('branches'). 

On the other hand, if we consider the two RH neutrino model as a limit of the 3 RH 
neutrino model for M3 ^ 10^^ GeV, then the orthogonal R matrix tends to 

(0 cos z C sin z \ 
-sinz Ccosz j (51) 

10 / 

and 

(cos z C sin 2; \ 
- sin 2 C cos ^ . (52) 

01/ 

Notice that the two branches cannot be obtained from each other with a continuous 
variation of the complex angle. This can be clearly seen if one considers the following 
general parametrization of the orthogonal matrix as a product of three complex rotations, 

-R(2;23, 2;i3, ^12) = ('-^23(2:23) -^13(2^13) -^12(^12) , (53) 

12 



1 











COS 223 


sin 223 





- sin Z23 


cos 223 



where 

' "" '^ n \ / (.Qg 213 sin 213 \ / cos 212 sin 212 \ / r A\ 

-R13 - 1 , i?-12 - -sin 212 cos 212 (54) 

y —sin 213 ^ cos 213 / V '^ ^ ^ J 

and where the overall sign (' = ±1 takes into account the possibility of a parity transfor- 
mation as well. Within this general case the two RH neutrino model R matrix for NH 
eq. ( I5TI) is obtained for 213 = z, 223 = C^i2 = 7r/2 and (' = (, clearly showing that the 
two branches for ( = ±1 cannot be obtained from each other with a continuous variation 
of the complex angle z (analogously for IH). 

We will refer in the following to this kind of view of the two RH neutrino model. We 
can perform scans over z for a fixed U, mi,Mi. 

Neutrino oscillation experiments measure two neutrino mass-squared differences. In 
the case of the two RH neutrino model for NH one has mi = 0, m2 = mgoi = v^"^^ = 
(0.00875 ± 0.00012) eV and m^ — ?Ti| = Aml^^. The heaviest neutrino has therefore a 



mass ms = matm = v ^"^atm + ^"^soi ~ (0.050 ± 0.001) eV [3l|- In the case of IH one 
has ma = 0, m2 = matm and mi = \/ml^^ - ml^^ 

We will adopt the following parametrization for the matrix U in terms of the mixing 



angles, the Dirac phase 5 and the Majorana phase 021 32| 

C12 Ci3 S12 Ci3 Si3 e ' ' 

U = -S12 C23 - C12 S23 si3 e"^ C12 C23 - S12 S23 Si3 e"^ S23 Ci3 

\ -312 523 — 012 023 313 6' — C12 S23 — S12 C23 S13 c' C23 C13 J 



and the following 2 a ranges for the three mixing angles 31 1 



■ diag (l, e' "2' , 1 j 
(55) 



^12 = (31.3° -36.3°), ^23 = (38.5° - 52.5°) , ^13 = (0° - 11.5°) . (56) 

As we will see, there will be some sensitivity to the low energy neutrino parameters, 
in particular to the value of ^13, of the Dirac phase and of the Majorana phase. We will 
therefore perform the scans with the following 4 benchmark Upmns choices A,B,C and 
D: 

A: ^13 = 0, 5 = 0, ^ = (57) 

B: 9r3 = 11.5°, 6 = 0,^ = (58) 

C: 9rs = 11.5°, 6 = ^,^ = (59) 
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D: 013 = 11.5°, 5 = 0,^ = ^, (60) 

where for all benchmarks the solar mixing angle and the atmospheric mixing angle are 
fixed to 6*12 = 34° and ^23 = 45° which are chosen to be close to their best fit values. Notice 
that benchmark A is close to tri-bimaximal (TB) mixing, with no low energy CP violation 
in the Dirac or Majorana sectors, while the remaining benchmarks all feature the highest 
allowed reactor angle consistent with the recent T2K electron appearance results 33|. 
On the other hand, varying the atmospheric and solar angles within their experimentally 
allowed ranges has little effect on the results, so all benchmarks have the fixed atmospheric 
and solar angles above. Benchmark B involves no CP violation in the low energy Dirac or 
Majorana sectors, with any CP violation arising from the high energy see-saw mechanism 
parametrized by the complex angle z. Benchmark C involves maximal low energy CP 
violation via the Dirac phase, corresponding to the oscillation phase 5 = 7r/2, but has 
zero low energy CP violation via the Majorana phase, with a2i/2 = 0. Benchmark D 
involves maximal CP violation from the Majorana sector, 021/2 = 7r/2, but zero CP 
violation in the Dirac sector, 6 = 0. These benchmark points are thus chosen to span 
the relevant parameter space and to illustrate the effect of the different sources of CP 
violation. 

4.2 Decay parameters in the orthogonal parametrization 

We can start first expressing the quantities rhi, Ki, Ki^ and pi2 in the orthogonal 
parametrization. For the effective neutrino masses and the total decay parameters one 
has simply 

fhi = y^ mk\Rik? , and i^i = V" — \Rik? ■ (61) 



m^, 



Substituting eq. ([78]) for m/j^, = y/Wi ^^ y/rn^U^k Rik into Ki^, = \ mu^, \ /{Mi m*) one 
obtains 

1 



K — 

m* 



7 ^ \/^>Thc Uak Rik 



(62) 



From this expression and from the definition of Ki^ in eq. (I13p . one then also obtains 



|2 



pO ^ \Y.k'V^^Uak'Rik'\ .gg. 
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Finally, substituting eq. f l75]) into eq. (EI]) for pi2 and simplifying via J2a ^ak ^ak' = ^kk' 
yields 

1 



Pl2 



nil ^^2 



(64) 



7 ^ "^fc -Rifc -R2fc 

With two RH neutrinos, we may express all quantities explicitly in terms of complex 
angle z for NH [mi = 0) as 

mi = msoi I cos^l^ + matm I sin2;|^, /fi = /Csoi | cos2;p + /fatm | sinz|^ , (65) 

"^2 = "^soi I sin^l^ + matm Icoszl^, K2 = Ksoi\smz\'^ + Ks.tmlcoszl'^ , (66) 

1 2 

Kia = —- ly/nisoi Ua2 COS z + ( y/m^t^ 11^3 slu z\ , (67) 

m* 

1 2 

i^2Q = — - ICVm^Uas COS 2; - y/rn~iUa2 smz\ , (68) 

m* 



and 



Pi2 = -I — ^ |"^atm COS 2: siu* z — mgoi COS* z sin 2; I , (69) 

mi m.2 

where we defined Ksoi = rnso\/m,i, ~ 10 and i^atm = '^atm/'Tiv, ~ 50. 

For IH, (ms = 0), we can approximate mi ~ m2 = m.atm and simplify further 

mi ~ m2 — matm (I cos zp + I sinzp) = matm cosh[2Im2;] , (70) 

i^i ~ i^2 =^ A'atm (I cos2;p + | sinz^) = K^tm cosh[2Im2;] , (71) 

Kia ^ -ft'atm \Ual COS Z + ( Ua2 Sm zf , (72) 

K2a - -ft'atm |C Ua2 COS Z - Ual Sill zf (73) 

and 

4 llm [cos* 2;sin2;l| , , .„^ ,,2 ,„,, 

^12 ^ ^ ,2 I • I 2 = |tanh[2Im2;]|^ . (74) 

(I coszf + I sm2;|"'j 

In Fig. [1] we show contour plots of the flavoured decay parameters i^i^,i^i^ and 
-^27; K2T in the relevant region of the z complex plane for NH and for the benchmark 
case B, since this will prove the case maximizing the effect of the N2 asymmetry pro- 
duction. Notice that Fig. [1] is periodic in vr along the Re^; axis as can be also inferred 
analytically from eqs. ( 1671) and ( l68l) using double angle identities. The most significant 
feature to be noticed at this stage is that for most of the parameter space Kia ^ 1 
holds. In these regions a strong wash-out regime is realized and this implies that the 
dependence of the results on the initial conditions is negligible and corrections due to the 
effects that we have listed earlier, after the kinetic equations, are at most 0{1) factors. 
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Figure 1: Contour plots showing Ki^ (upper left panel), Ki^ (upper right panel), 7^27 
(lower left panel) and K2T (lower right panel) dependence on complex angle z for bench- 
mark B (cf. eq. (EHD), C = +1, and NH. 

On the other hand, as we will discuss in section [5l there are two interesting new favoured 
regions for leptogenesis around z± ~ 7r/2 for NH, where the decay asymmetry from N2 
decays dominates over the one from A^^i decays (W2-dominated regions'). Fig. [1] shows 
that in this region Kia ~ 2 -^ 5. We are therefore in a 'optimal washout' region where 
thermal leptogenesis works most efficiently but still the dependence on the initial condi- 
tions amounts not more than ~ 50%. We have therefore decided to show the results just 
for the case of initial thermal A^2-abundance since the expressions for the efficiency factors 
are particularly simple. Very similar results are obtained for the other benchmark cases 
as well. 

At the same time, with Ki^ and Kir 3> 1, the asymmetries At- and A^^ produced 
from A^2-decays are efficiently washed out by A^^i inverse processes, and practically only 
the orthogonal component Ax, with size determined by 1 — Pi2, survives. Fig. |2] shows 
the contour plot of pu which indicates that the quantity significantly differs from unity 
in general. For NH pi2 is periodic in vr along the Kez axis and is approximately periodic 
in 7r/2 as can be inferred from eq. ( 169|) . Notice also that for IH, pi2 depends on Im^; only, 
in agreement with the analytical expression eq. (17^ . One can already see that in the new 
favoured regions around z ~ ±7r/2 the quantity 1 — pi2 is maximal. We thus find good 
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Figure 2: Contour plots of pi2 for NH (left panel) and IH (right panel). 



conditions for leptogenesis regarding washout from N2 as well as from Ni processes. 



4.3 CP Asymmetries in the orthogonal parametrization 

Let us now re-express the CP asymmetries in the orthogonal parametrization. The ex- 
pression ( IT6|) for the CP asymmetries can be recast as 



16 7rt;2 {rn}j^mo) 



r^Ef-; 



i¥^i 



iW/Mf) 



' >J i-i 



M,/M, ^'^ 3(M2/M2-1) 






(75) 
(76) 



in an obvious notation where we have defined, 



^ = Im 



[m^D)iai"^D)^Am^D"^D 



Jt, = Im 



{"lri^D)^o.^"^D)^j{m}^VflD)J, . (77) 



It is evident that X"- = — X°- and JZ" = —Jfi- In terms of the R-matrix we write from 
eq. (USD, 

mn = UD^,R^D^. (78) 

Then using this we find: 



I^j = Mi Mi 



M E 



.mkruk 



1) rnk'/U^kUak' R*kRjk'R*k" Rjk" 



k,k'k" 



j;« = M,M,Im[ Yl {rnurnu,f^mu„Ul^U^u'R*kR3k'R]k"R^k" 



k.k'.k" 



(79) 



^0) 



In order to simply the notation, it will prove convenient to introduce the ratios 



' in 



e{M, 



la X 

r 

1 ' %a 



e{M, 
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ia ^J 



-J 



e{M,) ■ 



^1) 



where 



£(Mi 



3 Ml matm 



IQlT v^ 



10" 



M, 



IQio GeV 



^2) 



is the upper bound for the total CP asymmetries 35|] that is therefore used as a reference 
value. 



4.3.1 Lightest RH neutrino CP asymmetries 

We can start from the lightest RH neutrino CP asymmetries e\a. We first write them 
including the third heaviest RH neutrino corresponding to the terms j = 3. We need then 
to specialize the general expressions above for Tjj and Jij to the case i = \ obtaining 



Xy = MiMjIm|^ ^ {mkmk'y/'^mk»U^f.Uak'R*ikRjk'R*ik"Rjk" 



k,k'k" 



and 



J^. = MiMjlm\ Y, {mkmk>Y'^mk"UlkUak' RlkR3k'R*k" Rik" 



k.k'.k" 



^3) 



(84) 



When we sum over j in the first term of the eq. fl75|) for i = 1 containing ly, thanks to 
R orthogonality and considering that for M2 > 3 Mi we can approximate ^{Mj/Mf) ^ 1. 
Then, only terms k' = k" survive and one obtains 11 1 



"la - 7 , 



k,k' 



fTik' y/ruk' ruk 

ffli matm 



\in.\Uak U^k' Rlk Rlk' 



^5) 



where the effective neutrino masses rrii can be written in terms of the R-matrix using 
eq. (1611) . This term is bounded by [lO| 

maxj [rrii] 



K < \Pl 



la 



'^atm 



6) 



and it is the only term that has been considered in all previous analyses of leptogenesis in 
the two RH neutrino model so far. It is useful to give a derivation of this upper bound. 
One can first write 



'la — 



< 



max,rra,- 



mi matm 



max, m,- 



mi matm 



^ y/mJJm^lm [Uak U*^' Rik Rik' 



k,k' 



^7) 



7 ^ y/mk Uak Rlk 



YVm^Kk'Rik' 



and then, using the eq. (!63|) and defining 



pO _ \T.k'V^^Uak'Rlk'\ , ^ 
-' la — ~ ^ -■-' 



mi 
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^9) 



one finally obtains 



r 



z I 

la 



niax,m,' 



<—^\IPLPL- (90) 

Notice that, if P°„ < 1, then P°„, and consequently rf„, can saturate its upper bound only 
if f/ is not real, i.e. only if the low energy phases are non-vanishing. This is because, if 
V is real, then one obtains P^^ = P^^ and the CP asymmetry would be necessarily below 
the upper bound. One can object that, if P^^ = 1, then the bound is saturated even if U 
is real. This is true for the CP asymmetry but it is not the case that maximizes the total 
asymmetry and that allows to saturate the lower bound on Mi. Let us explain briefly 
why. 

The possibility to relax considerably the lower bound in the usual flavoured A'^i dom- 
inated scenario compared to the unflavoured case relies on the so called one-flavor domi- 
nated scenario. In this scenario for a certain flavour a one has Pi^ <^ 1 while at the same 
time the upper bound on |rf^| is saturated. In this way the washout is reduced such that 
the efficiency factor gets enhanced approximately as oc 1/(P°q,)^'^ while the CP asymme- 
try is suppressed only oc \^Pi^ so that in the end the final asymmetry is approximately 
enhanced oc 1/(P°,)0-^. 

This derivation clearly shows, in general, that the possibility to have maximal asym- 
metry from the term rf^, necessarily relies on non-vanishing low energy phases in a way 
that PfL^ 1 ^ Pia- This analytical general demonstration is con&med by the numerical 



results [28 1 



For the other term the situation is quite different. The second term containing J'y 



cannot be simplified using the R orthogonality and one obtains 16| 



i»^~o Z^ TF — '^r~:^ lm[U„f^Uak' RikRjk' Rjk"Rik"\- (91) 

'-' ■ 1 1 1 1 II ^^■'■i III' 111'' atm 



r 

j,k,k',k" 



Let us now specialize the expressions eqs. ( 185|) and (lOTl) for r[^ and r^ to the two RH 
neutrino case using the special forms for the orthogonal matrix R (cf. (1511) and fl52p ) for 
NH and IH respectively. One can immediately check that the j = 3 terms vanish and for 



NH one obtains 11 1 

r?„ = '^lH^in'z\(\U^,?-^\U^A (92) 

mi \ mi^^ J 

+ C ~'°' ^ {("latm - "Isol) Im[f/a2 U^.^] Re[sin Z COS z] 

mi matm 

+ ("^atm + "^sol) Re[t/Q2 f^as] Im[sin Z COS z] } 
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and 

-^- = ^^{if'"t^^"'^niM^-if^-n (93) 

+ C —:::-^^ — — [(matm " "^soi) Im[f/*2 Uas] Re[sin z COS* z] (| COS z|^ + | sin ;z| 

?Til matm 

+ (''Tiatm + ■'Tisoi) Re[f/*2 f/os] Im[sin 2; COS* z] (I COS 2;p - | sin 2;^)] } . 
In terms of Kez, Iraz the dominant term rf^ is 



2^ 



'^atm 1 . rr^T-. n • i rr^r i I \tt i2 '^i 



mi 2 \ m: 



rf^ = -^ - sin[2Rez] sinh[2Imz] \U^^\' - -^ |f/,2r (94) 



2 

''atm 



+ - C ^^^-^ — ^ {("^atm - "^soi) Im[t/a2 ^^s] sin[2Imz] cosh[2Im2;] 
2 mi matm 

+ ("^atm + "^soi) Re[[/„2 t/*3] cos[2Re2;] sinh[2Imz]} . 
Analogously, for the case of IH and approximating mi ^ m2 — ^^a.tn\i one obtains 

T '^atm 



1" mi 



{Imfsin^^] (|f/„ip - |f/„2n -2CRe[t/„it/*2]Im[sin^ cos^]} (95) 



and 



< - ^^^ Wsin^^](|t/„ip-|t/„,n (96) 

+ 2C (I sinzp - I cos2;p) Re[t/Qi t/*2] Im[sin2 cos*^]} . 

Notice that while the terms rjj^ are proportional to M1/M2, the terms r^^, are not. In 
Figure [3] we have plotted the quantities rf„ and r'^a/^'ia = ^la/^'ia fo^ ^^e benchmark B 
UpMNS choice eq. f l58|) . C = +1 and M2/M1 = 3. Once again there is periodicity in tt 
along Re^;, for the same reasons as with Figs. [T|2l One can notice how r^al^fa ^ ^ 
both for a = 7 and a = t. Only in a very fine tuned region this ratio gets up to about 
0.5. Therefore, it will prove out the term rj^ can be safely neglected in the regions of 
interest for this study. Also, one can notice that r^^,, the dominant contribution to the 
baryon asymmetry from A'^i decays, is suppressed in the region z ~ -7r/2, hence this region 
is potentially dominated by A'^2 decays. 

From the lower panels of fig. |3] one can notice how T^al'^\a i^ independent of Im[2], 
eq. f l58|) . This is because Im[f/Q,2 f/*3] = such that the middle terms vanish from eqs. fl92l) 
and fl93|) . It can then be shown that all the dependance of r\^ and r^ on Im[2;] is contained 
in the common factor sinh[2Im2;] jm\ which then cancels in the ratio Ti^/r^^. 
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Figure 3: Contour plots of r\^ (upper left panel) , r[^ (upper right panel), \T^-~^lr[^\ = 
\^i"//^%\ (lower left panel) and |r^/rj^| = |5j^/ef^| (lower right panel) for NH, benchmark 
B eq. (EZD, C = +1 and M2/M1 = 3. 
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Figure 4: Contour plots of r^ (upper left panel), r^ (upper right panel), |?"^/Tf^| = 
|e^/ef^| (lower left panel) and |?"2^/rf^| = k^/^^frl (lower right panel) for benchmark A, 
eq. (ET]), C = +1, M2/M1 = 3 and NH. 
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4.3.2 Next-to-lightest RH neutrino 

Let us now turn to consider the case i = 2. For j = 3 we have 

1^3 = M2 Ms ^ TUk" ^Jmkmk>lm[UlfJJak' R*2k -Rsfc' Rlk" Rsk"] 



k,k'k" 



and 



J23 = M2M3 ^ TJik" ^rriki ruk lm[U^i^Uak' Rlk R-sk' Rlk" R^k" 



(97) 



(9^ 



k,k'k" 



It is easy to check that both two terms vanish in the two RH neutrino case. 
On the other hand the two terms for j = 1, 

X21 = M2M1 ^ rrik" y/mkmk'lm[U^f^Uak' Rlk Rik' Rlk" Rik"] , 

k,k',k" 



and 



J21 = M2M1 ^ mfc" i/mfc/ mk lm[UlfJJak' Rlk Rik' Rlk" ^2fc"] 



(99) 



(100) 



k,k'k" 



do not vanish and they lead, in the hierarchical limit M2 > 3 Mi, to final values for r 



2a 



and rgQ, given respectively by 



T 
'2a — 



4 /Ml 



3 VM2 






and 



7 

'2a — 



E 



k,k'k" 



nik" ^rrik' ruk 



rrik" y/mk my 



m2 matm 



lm[f/*^.t/Qfc/ i?2fc Rlk' Rlk" Rlk"] , 



3 ,^„ 1^2 matm 

k,k',k" 



lm[[/*fc f/ofc' -R2A: -^Ifc' -Rlfc" -R2fc"] • 



(101) 
(102) 



This second term r^ will clearly tend to dominate over rf^, oc (M1/M2). However, since 
the dependence on the complex parameter z is different, one cannot exclude that, in 
some region of the parameter space, r^^, can give a non negligible contribution. We have 
therefore safely taken into account this term checking indeed that is negligible. 
If we specialize the expressions to the two RH neutrino model we obtain for NH 



2a 



4 Ml 
'3 M2 



'■>'t'- 



m, 



atm 



mi 



Im[sin^ z] 



\U^ '^ 



m,' 



^aSl 



sol \TT 1 2 



m 



'o2\ 



atm 



(103) 



\/ i\\XX\ sol / \r -j^nr j-nri lO 1 iOt 

+ C ~ ("^atm - "T-soi) \m\U ^2 f^sJ Re[sin 2; COS zj [I COS z| + I sin 2; I J 



+ C 



m2 matm 

A/m.atm m-sol 

m2 matm 



'"^atm + '^soi) Re[t^a2 f^os] Ini[sin z COS* z\ [I COS zp — I sin z I 
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and 



r-^ _ ^ "^'^"ii^rci^s^i rirr J2 _ irr j2 



lm[sm^z] [\Ua2\^ - \Ua3\^] (104) 



3 7712 

+ 77 C ~ ''^"' — ^ ("^atm - "^soi) Im[f/*2 ^os] Re[sin z COS* z] [| COS zp + I sin z\'^] 

3 7712 matm 

\/ flytm sol / \ rj- ^ r -j-nri lOi iOt 

+ o'^^ (matm + "^soi)Re[f/„2f43jIm[sinz COS zj [I cos2;| -|sinz|-J 

3 ?Tl2 ?Tlatm 

In terms of Re^, liaz the dominant term r^ is 



2q 



3 m2 

3 7712 matm 



sinpRe^] sinh[2Im^] [|?7a2|^ - It/'asP] (105) 



C ~ ''^'" — ^ ("^atm - "^soi) Im[f/*2 f^s] sin[2Re2] cosh[2Imz] 



- C ~ ^^"^ — — ("^atm + "^soi) Re[t/*2 ^as] cos[2Re^] sinh[2Imz] 



"^2 "^atm 



For IH we obtain 



"^atm rT_r„;„2^i /\tt |2 \tt |2 



2a 



{Im[sin2^] {\U^i\^ - \Ua2\^) (106) 

mi 

2C (I sin z\^ Im[U^-^ Ua2] Re[sin2; cos* 2] 
I cos z\'^ Re[U*i Ua2] Im[sin z cos* z]) 



and 



^J _ ^atm rT_r.;^2,i Htt |2 Irr |2 
2a ~ 



{Hsin^^] (|?7„i|2-|?7„2r) (107) 

7712 

+ 2 C Re[f/Q,i [7*2] Im[sin 2 cos* 2;] [| sin z\'^ - | cos z\'^] } . 

In Figure m we have plotted r^^ and r^/r^^ for C = +1, M2/M1 = 3, benchmark Upmns 
choice A (c.f. eq. (15 7|) ) and NH. This time, as one can see from the figures, one has 
{^2al^2a) ^ ^ ^^^ ^^^ valucs of z and M1/M2 (since '"^/^fo, gets even larger if M2/M1 > 3 
), implying that the term r^ dominates and that r^^, can be safely neglected. It can again 
be seen in fig. Hlthat ri^^/r^a depends only on Res; for the same reasons as with T^Jr\^. 
Once again there is periodicity in vr along Res;, for the same reasons as with Figs. [T|2|3l 

Crucially, we find that rg^, the dominant contribution to the baryon asymmetry from 
N2 decays, is maximised in the regions z ~ ±7r/2 (just above and below the Ims; = 
line), in contrast to rj^,, the dominant contribution from A'^i decays, which is minimised in 
this region (see fig [3]) . Given this result and the favourable values of 7^27 and pi2 , shown 
in fig [T] and fig [2] respectively, one expect the z ~ ±7r/2 regions will be N2 dominated. 

Notice that we have not shown any figure for the case of IH since it will turn out that 
the contribution from the next-to-lightest RH neutrinos to the final asymmetry is always 
negligible. 
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5 Constraints on the parameter space and A^i versus 
N2 contribution 

We can now finally go back to the expression for the final asymmetry (cf. eqs.f l40p . (141 p 
and fH2|) ) and recast it within the orthogonal parametrization. This can be written as the 
sum of four terms, 



A^^ 



B~L 



A^,^) 



A^,J) 



rf(2,2:) 



N'^^l>^m^i> + N'^^r' + N- 



rf(i) 



f(2,J) 



:io8i 



The sum of the first two terms is the contribution N^\ from the lightest RH neutrinos, 
<i'J)(z,f/,Mi) = £(Mi) Y, rUz,U)K{K,^), (109) 

iVjSi'f(z,f/,Mi,Mi/M2) = £(Mi) J2 rt{z,U,MjM2)K{K^o.), (HO) 



a=r,7 



and it should be noticed that only the second one depends on M2. 

Analogously the sum of the last two terms in eq. (llUSp is the contribution N-^J^ from 
the next-to-lightest RH neutrinos, 



A2,x) 



N'^tl\z,U,M^,M^/M2) 



N 



B-L > 



^,t/,Mi) 



e{M^)f^{z,U), 



f{z,U) 



where 



Ml 
M2 



In 



M2 

Ml 



f{z,U) = pi2rj^K,iK2^)e~ 



-Ki^ 



(111) 
(112) 

(113) 



+ {1 - P12) rL k{K2^) + r^^ K{K2r) e- 



-Klr 



and 



f^{z, U) = P12 rl t,{K2,) e-'-f""^- + (1 - p^2) r^ <K2,) + ri /.(Ka.) e""^^- . (114) 

Notice that this time the first term depends on M2 while the second does not. We can 
then write the total asymmetry as 



A^Ll 



[A^2-? + Nf_:i^]{z, U, Ml) [1 + ^^1 + 62]{z, U, Ml, M1/M2 



(115) 



f(i.^)/r/\rf(i.2:) 



B-L 



N 



m,Jh 



.f(2,J) ,r,.f(l,J) 



^_,,and62^Nr:t>/[N'ti' + N's.L 



Tm,Jh 



We 



where we defined 61 = N^^^'/[N^ 

found that ^1,^2 < 0.05 for any choice of M1/M2, -z, ?7. Therefore, one can conclude 
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that the total final asymmetry is independent of M2 with very good accuracy. It should 
be however remembered that our calculation of N^\ holds for M2 < 10^^ GeV and 
M2/M1 > 3, implying Mi < (100/3) x 10^° GeV. As such, when N2 decays are included 
the largest value of Mi we will allow is Mi = 30 x lO^^GeV, whereas when N2 decays are 
neglected, we will consider values as large as Mi = 100 x 10^°GeV. 

In Fig. [3] we show the contours plots for Mi obtained imposing successful leptogenesis, 
i.e. 7] = 7]'^^^ (we used the 2a lower value rf^^^^ = 5.9 x 10"^*^), for ( = +1 and for initial 
thermal abundance. The four panels correspond to the four benchmark cases A, B, C 
and D in the NH case. The solid lines are obtained including the contribution N^_^ in 
the final asymmetry and therefore represent the main result of the paper. These have to 
be compared with the dashed lines where this contribution is neglected. In Fig. [5] and 
indeed in all subsequent figures, one can notice again a periodicity in tt along Kez. This 
is because final asymmetries are given from eq. (llOSp . for which all terms are dependant 
upon quantities periodic in vr along Re^ (these quantities being the washouts, pu and the 
CP asymmetries). As one can see, on most of the regions leptogenesis is A'^i-dominated 
as one would expect |j. However, there are two regions, around z ~ ±7r/2, where the 
asymmetry is A''2-dominated. If N^_^ is neglected, this region would be only partially 
accessible and in any case only for quite large values Mi > 30 x 10^° GeVQ 

When the contribution N^_^ is taken into account one can have successful leptogenesis 
for Ml values as low as 7 x 10^° GeV for benchmark case B and initial thermal A^2- 
abundance. The existence of these W2-dominated regions' is the result of a combination 
of different effects: i) the value of (1 — pu), setting the size of the contribution from A'^2 
decays that survives the A'^i washout, is maximal in these regions as one can see from Fig. [21 
ii) the wash-out at the production is in these region minimum as one can see from the plots 
of K2r and -ft'27 (cf. Fig. [1]); iii) the A'^2-flavoured CP asymmetries are not suppressed in 
these regions contrarily to the A'^i flavoured CP asymmetries. It is interesting to compare 
the results obtained for the 4 different benchmark cases. A comparison between the case 
A (upper left panel) and the case B (upper right panel) shows that large values of ^13 and 



■^The A'^i-doniinated regions are approximately invariant for z — > — z, implying N^^j^{z) ~ 
N^^j^ {—z). This is because rf^ is dominated by the first term in the eq. (|M)) . exactly invariant for 
z — > —2, and because the Kia are also approximately invariant for z — > —2 (see upper panels in Fig. 1). 

"* Notice that this time there is no invariance with respect to 2 — > —2 since 7*2^ is dominated either 
by the third term (for cases A, B, C,) or by the second term (for case D) in eq. (|105p that are not 
invariant for z ^ —z. On the other hand the second term is invariant for Re 2 -^ —Re 2 and therefore 
one could naively expect a specular region at 2 ~ —tt/2. However, notice that A'27 is not invariant for 
Re 2 — > —Re 2. In this way, for negative Re 2 and same values of \z\, the wash-out is strong and prevents 
the existence of this specular region. 
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Figure 5: Contours plots in the 2;-plane of the Mi values obtained imposing successful 
leptogenesis (jib = Vb^'^) fo^ the NH case, C = +1 and benchmarks A (top left), B (top 
right), C (bottom left) and D (bottom right) fixing U. The solid lines are obtained taking 
into account the contribution N^_^ to the final asymmetry while the dashed lines are 
obtained neglecting this contribution. Contours are labelled with the value of Mi in units 
of 10i°GeV. 
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no Dirac phase enhance N^_^ so that the A^2-doniinated regions get enlarged. On the 
other hand a comparison between B and C shows that a Dirac phase seems to suppress 

1(2'] 

N^_^. A comparison between B and D shows that a Majorana phase seems just to 
change the position of the A'^2-dominated regions without consistently modify their size 
differently from the A^i-dominated regions that are instead maximized by the presence of 



non- vanishing Majorana phase as known 13|, |28|. Interestingly, for non-zero Majorana 



phase the new region where leptogenesis is favoured now overlaps with the Im(2;) = 
axis. This means that CP violation for A'^2-dominated leptogenesis can be successfully 
induced just by the Majorana phase. We have also checked that varying the low energy 
parameters within the ranges of values set by the 4 benchmark cases, one has a continuous 
variation of the allowed regions. 

It can be seen that the A'^2-dominated regions are maximal in case B. For this reason 
in Fig. [6] we show a zoom of the iV2-dominated regions around z = 7r/2 for case B. This 
figure represents one of the main results of this paper. Notice that if we were starting 
with initial vanishing abundance the A^2-dominated regions just slightly shrink and the 
lower bound on Mi becomes just slightly more restrictive, 10 x 10^*^ GeV instead than 
7 X 10^° GeV. 

On the other hand if we consider the IH case, the situation is very different as one 
can see from Fig. [71 The much stronger wash-out acting both on the A'^i and on the A'^2 
contributions suppresses the final asymmetry in a way that large fraction of the allowed 
regions disappear, including the iV2-dominated regions. The surviving allowed regions are 
therefore strongly reduced and strictly A'^i-dominated. Analogous results are obtained for 
the branch ^ = — 1, shown in figures [H A comparison between the plots obtained for the 
two branches shows that the the finally asymmetry is invariant for (^, z) — ;■ (— ^, —z) and 
this is confirmed by the analytical expressions both for the flavoured decay parameters 
determining the wash-out and for the CP asymmetries. 

6 Leptogenesis from two RH neutrinos in models with 
Light Sequential Dominance 

In section [5] (c.f. Figure E]) we have seen that two new favoured region for leptogenesis 
have appeared where z ~ ±'7r/2, for (^ = ±1, and for NH. Compared to previous studies 
where the production of the baryon asymmetry in this region of parameters was thought 
to be very suppressed, we found that, due to effects from A^2 decays, leptogenesis is quite 
efficient and can be realised with comparatively low Mi of about 7 x 10^*^ GeV. This 
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Figure 6: Contours plots in the 2;-plane of the Mi values obtained imposing successful 
leptogenesis (77^ = Vb^^^) fo^ the NH case, C = +1- Enlargement of the benchmark B 
case from previous figure. This is a particularly interesting case, since it maximises the 
N2 dominated region around z ~ zlsd = 7r/2. 
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Figure 7: Contours plots in the 2;-plane of the Mi values obtained imposing successful 
leptogenesis (77^ = rj^^^^) for the IH case, C = +1 and benchmarks A (top left), B (top 
right), C (bottom left) and D (bottom right) fixing U. The solid lines are obtained taking 
into account the contribution N^_j^ to the final asymmetry while the dashed lines are 
obtained neglecting this contribution. Contours are labelled with the value of Mi in units 
of lO^^GeV. 
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Figure 8: Contours plots in the 2;-plane of the Mi values obtained imposing successful 
leptogenesis (77^ = r/^*^^) for the NH case, ( = —1 and benchmarks A (top left), B (top 
right), C (bottom left) and D (bottom right) fixing U. The solid lines are obtained taking 
into account the contribution N^_j^ to the final asymmetry while the dashed lines are 
obtained neglecting this contribution. Contours are labelled with the value of Mi in units 
of lO^^GeV. 
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Figure 9: Contours plots in the 2;-plane of the Mi values obtained imposing successful 
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Figure 10: Contours plots in the ^-plane of the Mi values obtained imposing successful 
leptogenesis {jib = Vb^^^) ^^t the IH case, ( = —1 and benchmarks A (top left), B (top 
right), C (bottom left) and D (bottom right) fixing U. The solid lines are obtained taking 
into account the contribution N^_j^ to the final asymmetry while the dashed lines are 
obtained neglecting this contribution. Contours are labelled with the value of Mi in units 
of 10i°GeV. 
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result is particularly interesting since z ~ ±7r/2 corresponds to the class of neutrino 
mass models with Light Sequential Dominance (LSD) [2|, as we will discuss below. The 
dictionary between the parameter z and the Sequential Dominance (SD) parameters will 
be given explicitly in section 16.21 Finally, in 16.31 we will discuss the decay asymmetries in 
an explicit example scenario of LSD and show the enhancement of the asymmetry from 
A^2 decays analytically in terms of SD parameters and the deviation from TB mixing. 



6.1 Light Sequential Dominance 

In models with SD, the RH neutrinos contribute to the neutrino mass matrix with "se- 
quential" strength, leading to a NH. In LSD, the lightest RH neutrino A^^i provides the 
largest "dominant" contribution, whereas the second lightest RH neutrino contributes 
subdominantly. When the heaviest RH neutrino (almost) decouples, we arrive (approxi- 
mately) at a 2 RH neutrino model. 

To understand how SD, and in particular LSD works, we begin by writing the RH 
neutrino Major ana mass matrix Mrr in a diagonal basis as 

/ Ml \ 
Mrr= M2 , (116) 

V Mg / 

where M\ < M2 < M3. In this basis we write the neutrino (Dirac) Yukawa matrix X,y in 
terms of (1,3) column vectors Ai, Bi, Ci as 



K 



ABC 



;ii7) 



in the convention where the Yukawa matrix is given in left-right convention. Explicitly 

we have 

/ Ai Bi Ci \ 

Y,= A2 B2 C2 . (118) 

V A3 53 C3 / 

The Dirac neutrino mass matrix is then given by m^^^^ = X^v^. The term for the 
light neutrino masses in the effective Lagrangian (after electroweak symmetry breaking), 
resulting from integrating out the massive RH neutrinos, now reads 



c 



eff 



:ii9) 



Ml M2 M3 

where z/j {i = 1,2,3) are the left-handed neutrino fields. As stated above, LSD then 
corresponds to M3 — > 00 so that the third term becomes negligible, with the second term 
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sub dominant and the first term dominant [2l: 



Ma Mb Mc 

In addition, we shall assume that small ^13 and almost maximal ^23 require that 

l^il < 1^2! ~ \A2\. 
Constrained Sequential Dominance (CSD) is defined as 



Ai 


= 0, 


A2 


= 1^43 


Bi 


= \B2 


\^B 


= 0. 



IS, 



(120) 



:i2r 



(122) 
(123) 
(124) 
(125) 



CSD implies TB mixing 36| and vanishing leptogenesis if Mq » M^, Mb 12|, |37 |. 



6.2 An R matrix dictionary for LSD 

According to LSD, the "dominant" A'^i, i.e. its mass and Yukawa couplings, governs the 
largest light neutrino mass m^, whereas the "subdominant" A'^2 governs the lighter neutrino 
mass 777,2, while the decoupled A^3 is associated with mi — )■ 0. From eq. ( H71) it is then 
clear that, ignoring 7712/1^3 corrections, the R-matrix for LSD takes the approximate form 



/ 1 \ 



R 



LSD 



diag(±l,±l,l) 



(126) 



1 

VI 0/ 

where the four different combinations of the signs correspond physically to the four differ- 
ent combinations of signs of the Dirac matrix columns associated with the lightest two RH 
neutrinos of mass Mi and M2. The sign of the third column associated with M3 — )■ 00 
is irrelevant and has been dropped since it would in any case just redefine the overall 
sign of the Dirac mass matrix. These choices of signs are of course irrelevant for the 
light neutrino phenomenology, since the effect of the orthogonal R matrix cancels in the 
see-saw mechanism (by definition). The four choices of sign are also irrelevant for type 
I leptogenesis, since each column enters quadratically in both the asymmetry and the 
washout formulas in Eqs. ( TT6|) and ( TT3|) . independently of fiavour or whether A'^i or N2 is 
contributing. Comparing eq. (I126P to the parameterisation of R^^^^ for the 2 RH neutrino 
models in eq. (15T|) . we see that LSD just corresponds to 2; ~ ±7r/2 which correspond to 
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the new regions opened up by N2 leptogenesis that were observed numerically in the pre- 
vious section. To be precise the dictionary for the sign choices in eq. f ll26p are as follows: 
for the C = 1 branch, z ~ 7r/2, corresponds to diag(l, —1, 1), while z ^ — 7r/2, corre- 
sponds to diag(— 1, 1, 1); for the (^ = — 1 branch, z ^ 7r/2, corresponds to diag(— 1, —1, 1), 
while z K. — 7r/2, corresponds to diag(l,l,l). According to the above observation, all 
four of these regions will contribute identically to leptogenesis, as observed earlier in the 
numerical and analytical results (i.e. giving identical results for C, = ±1 and z ^ ±7r/2). 
We may expand eq. flSTl) for LSD for any one of these identical regions to leading order 
in m2/m^. For example consider the case C = — 1 and z ^ — 7r/2 corresponding to the case 
where all the Dirac columns have the same relative sign, diag(l,l,l). Then expanding 
eq. (15T!) around z ~ — vr/2, defining A f« 2 + |, we may write. 



R 



LSD 



/ A 


1 \ 


1 


-A 


[1 


/ 



(127) 



Using the results in [39[ we find useful analytic expressions which relate the R-matrix 
angle to the Yukawa matrix elements near the CSD limit of LSD corresponding to small 

A, 



Re(A) 



Re{A^B)v'^ 



1/2 .A/2 



{m3-m2)M^''M^ 



lm{A^B) 



2 



(m3 + m2)M^ M2 

Notice that A — )■ when A'^B — )■ to all orders in 7712/771^. This is just the case in 
CSD due to eq. ( 1125^ . Thus in the CSD limit of LSD eq. ( 1126P becomes exact |38| to all 
orders in 7n2/7n3. Clearly, leptogenesis vanishes in CSD which can be understood from 



the fact that the R-matrix in CSD is real and diagonal (up to a permutation) 37| or 



from the fact that A is orthogonal to B 12|. However in the next section we consider a 



perturbation of CSD, allowing leptogenesis but preserving TB mixing. 



6.3 Example: perturbing the CSD limit of LSD 



Using eq. flT6|) . we obtain, making the usual hierarchical RH neutrino mass assumption 
the A^^i contribution to the leptogenesis asymmetry parameter is given by: 
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Clearly the asymmetry vanishes in the case of CSD due to eq. (11251] . In this subsection 
we consider an example which violates CSD, but maintains TB mixing and stays close to 
LSD. 

Before we turn to an explicit example, let us state the expectation for the size of the 
decay asymmetries. We expect that, typically, 

3 maMi A 

^W^-J^^;^' ^le -^ ^^i/.,r. (130) 

The A'^2 contribution to the leptogenesis asymmetry parameter is given by the interference 
with the lighter RH neutrino in the loop via the second term in eq. (IT6l) . which is indeed 
often ignored in the literature: 

82. ^ -tI-t^Iiii [Bl{A^B)A^] . (131) 



This leads to typically. 



1671 B^B 



1 m^Mi Ai 



which should be compared to eq. (I130p . The A'^2 contribution to the decay asymmetries 
looks larger than the A'^i contribution. 

To compare the two asymmetries and the produced baryon asymmetry explicitly, let 
us now calculate the final asymmetries in a specific perturbation of the Light CSD form. 
As an example, we may consider 

(^1,^2,^3) = (0,a,a) (133) 

{BuB2,Bs) = {b,b + q,-b + q) (134) 

such that 

/ b Ci \ 

Y,= a b + q C2 . (135) 

\a -b + q Cs J 

Providing |g| ^ |6|, this perturbation of CSD but stays close to LSD and allows non-zero 
leptogenesis. Interestingly this perturbation of CSD also preserves TB mixing as discussed 
in |39|, where more details can be found. Note that z is given by eq. (I128p and therefore 
depends on a, b and q. 

For our example, we now obtain (assuming real a and neglecting q in B'^B): 

3 m2Mi\m\qb + q^] 3 msMi Im[-gfe + g^] 

^'^"^"16^^;^ Wb ' ^^^"^'m^^p Wb • ^^^^^ 
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The N2 contribution to the leptogenesis asymmetry parameter is given by the interference 
with the hghter RH neutrino in the loop via the second term in eq. ( lT6l) : 

2 m3MiIm[g6*] 



'^ 167r v^ B^B 

For the washout parameters, we obtain: 



e2r ~ -£27 (137) 



m 
and 



Ki^ = Kir^— (138) 



m 
The parameter pi2 is given by (neglecting q in the last step) 



K27 ~ i^2r ~ — • (139) 



^- " - AtA + BW " 6 ^''°^ 



For the final asymmetries from A^^i decay this means 

3 m2MiIm[g2 



whereas 



Nn « (.,, H- e,.HK,,) « -a^^^Kl™,/™-) (141) 

A'S « (1 - P,.K.K(/f..) « -5 J_!!!iMl!!^,(„,/,„.) , (142) 

So we can estimate: 

f(2) 

Nb-l _ m3 ft:(m2/m*) 5 \m[qh*] 
jy/(i) 77J2 K{m^/m*) 18 Im[g^] 

We see that, as already anticipated in the beginning of this subsection, there is an enhance- 
ment of the asymmetry from the N^ decay by a factor of ^^ (from the decay asymmetries) . 
Furthermore, there is another enhancement factor from the efficiency factor k given by 
K{m2/rn ) -Q^^Yi terms imply an enhancement of a factor of 5. On the other hand, there 
is a numerical factor of ?^ ^ which somewhat reduces the enhancement. Finally, the 
factor ^? 2} can get very large in the limit that q gets small, i.e. close to the CSD case. 
However, of course, closer to the CSD case the decay asymmetries get more and more 
suppressed. 

In summary, in models with Light Sequential Dominance (LSD) the asymmetry from 
the N2 decays is generically larger than the asymmetry from A^^i decays, in agreement 
with the results obtained in the previous sections in the R matrix parameterisation. We 
like to emphasise that in order to calculate the prospects for leptogenesis in models with 
LSD (in the two fiavour regime), it is thus crucial to include the A'^2 decays (which have 
previously been neglected). 
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7 Conclusions 

We have revisited leptogenesis in the minimal non-supersymmetric type I see-saw mech- 
anism with two hierarchical RH neutrinos (M2 > 3 Mi), including flavour effects and 
allowing both RH neutrinos Ni and N2 to contribute, rather than just the lightest RH 
neutrino Ni that has hitherto been considered. 

We emphasise two crucial ingredients of our analysis: i) the flavoured CP asymmetries 



have 



Deen calculated taking into account also terms that cancel in the total CP asymme- 



tries iJ] and that have been so far neglected within the two RH neutrino model; ii) Part 
of the asymmetry produced from A'^2 decays, that one orthogonal in flavour space to the 
lepton flavour ii produced by A''i decays, escapes the A''! washout |29|. 

Defining four benchmark points corresponding to a range of PMNS parameters, we 
have performed scans over the single complex angle z of the orthogonal matrix R for each 
of the two physically distinct branches C = ±1- For the case of a normal mass hierarchy, 
for each benchmark point we found that in regions around z ~ ±7r/2, the N2 contribution 
can dominate the contribution to leptogenesis. For benchmark B corresponding to a 
large reactor angle and zero low energy CP violation we found that the lightest RH 
neutrino mass may be decreased by about an order of magnitude in these regions, down to 
Ml ~ 7x 10^° GeV for initial thermal A''2-abundance, with the numerical results supported 
by analytic estimates. Other benchmarks with smaller reactor angle and/or low energy 
CP violating phases switched on exhibit similar results. 

It should be recalled that our results have been obtained in the strict limit M3 — )■ 00. 
Effects from a large but finite M3 could possibly enhance the impact of the N2 production 
enlarging the A'^2-dominated regions. 

These 7V2-dominated regions around z ~ ±7r/2 are quite interesting since they cor- 
respond to light sequential dominance in the hierarchical limit where the atmospheric 
neutrino mass m^ arises dominantly from the lightest RH neutrino of mass Mi, the solar 
neutrino mass 1712 arises dominantly from the second lightest RH neutrino of mass M2, 
and the lightest neutrino mass of mi is negligible due to a very large RH neutrino of mass 
M3. Such a scenario commonly arises in unified models based on a natural application of 
the see-saw mechanism 12|, |40|, |4l|, |42[ so the new results in this paper may be relevant 



to unified model building in large classes of models involving a NH. 
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